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Context. Recent studies have argued that galaxy mergers are not important drivers for the evolution of SO's, on the basis that mergers 
cannot preserve the coupling between the bulge and disk scale-lengths observed in these galaxies and the lack of correlation of their 
ratio with the SO Hubble type. However, about 70% of present-day SO's reside in groups, an environment where mergers and tidal 
interactions dominate galaxy evolution. 

Aims. We investigate whether the remnants resulting from collision-less N-body simulations of intermediate and minor mergers onto 
✓ ■ SO galaxies evolve fulfilling global structural relations observed between the bulges and disks of these galaxies. 

Methods. Different initial bulge-to-disk ratios of the primary SO have been considered, as well as different satellite densities, mass 
ratios, and orbits of the encounter. We have analysed the final morphology of the remnants in images simulating the typical observing 
conditions of SO surveys. We derive bulge+disk decompositions of the final remnants to compare their global bulge-to-disk structure 
with observations. 

Results. We show that all remnants present undisturbed SO morphologies according to the prescriptions of specialized surveys. The 
^ \ dry intermediate and minor mergers induce noticeable bulge growth (SOc — >S0b and SOb — >S0a), but affect negligibly to the bulge 

and disk scale-lengths. Therefore, if a coupling between these two components exists prior to the merger, the encounter does not break 
this coupling. This fact provides a simple explanation for the lack of correlation between the ratio of bulge and disk scale-lengths and 
the SO Hubble type reported by observations. 

Conclusions. These models prove that dry intermediate and minor mergers can induce global structural evolution within the se- 
quence of SO Hubble types compatible with observations, meaning that these processes should not be discarded from the evolutionary 
scenarios of SO's just on the basis of the strong coupling observed between the bulge and disk scale-lengths in these galaxies. 
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£NJ 1 . Introduction cently, locating SO galaxies in a sequence (SOa-SOb-SOc) paral- 
t— I ■ lei to the one of spirals (Sa-Sb-Sc), instead of placing them at 
Lenticular galaxies (SO's) are disk galaxies with a smooth ap- an intermediate location between ellipticals a nd spirals, as done 
pearan ce due to the i r lack of spiral arms and star-formin g re- in the traditional Hubble diagram (see L1 0; LaurikamenetalJ 
^ ! gions (iHubbldl 19261 Ide Vaucouleursll 19591: ISandagdfl96ll) . The 120111 : ICappellariet al kolU iKormendv & Bendet1l2012lh 
H. traditional picture of these galaxies as a transition class between The large dispersion of B/T for SO galaxies may be re- 
spirals and ellipticals is drastically changing in the last years, fleeting di fferent formation and evolution mechanisms, as sug- 
because although SO's tend to exhibit massive central bulges, gested by Ivan den Berghl dl976l) . The later appearance of SO's 
recent observations show that their bulge-to-dis k ratios (B/D) into the cosmic scenario and the noticeable rise of their popula- 
can sp an a wide range of values. In this sense, IBalcells et alj tion in clusters at the expense of spirals since z ~ 0.5 evidence 
(1200711 found no correlation between B/D and the typical se- that SO's must derive from them (|de Jondfl99"6l: ICourteau et al.l 



quence of Hubble types (SO-Sa-Sb-Sc): their SO galaxies had a 1996; Fa sano et al] 120001: iGrahaml 120 01b: Smith et al. 2005 



r 



mean B/D value similarto that of late-type spirals (B/D ~ 0.25). Postman et al. 2005; Balcells et alj|200 7; Poggi anti et alj|200i 

Studies using larger samples of galaxies including lenticulars 2009: IWilman et al.1 2009. L10). However, the main mechanisms 

show that, although the bulge-to-total luminosity ratio (B/T) in- after their formation and evolution are still unsettled, with data 

creases from late-type to early-type galaxies, some SO galaxies indicating that quite different processes may have conttibuted 

can exhibit B/T values similar to later types ( B/T « 0.1-0.3, see to it depending on the environment and on their luminosity 

IWeinzirletatll2009t lLaurikainen et all 120101 L10 henceforth). (Bar wav et alJl2007l 120091: iKannappan et alJl2009l:ISimard et al.l 

This fact has made three different groups independently to up- 2009; Huertas-Companv et al.ll2.01QI) . Thi s has bee n point ed out 

date the van den Bergh's version of the Hubble tuning fork re- in a recent review about SO galaxies by lAguerril d2012l) . who 
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derives this conclusion attending to the different properties of 
the photometric components of lenticular galaxies and their im- 
plications in the formation and evolution mechanisms of these 
galaxies. 

SO's, like most disk galaxies, show a strong coupling 
between the b ulge effective radius (r e ) and the disk scale 
length (h D , see lHunt et al.ll2004t lLaurikainen et~aT]l2009l L10). 
As argued by lLaurikainen et all the hierarchical merger sce- 
nario cannot account for such coupling, because bulges grow 
by mergers whereas disks grow later from the accretion 
of left-over gas, these two processes being thus decoupled 
(lLaurikainen et al.l 120091 L10). This reasoning has left merg- 
ers in the backstage of the SO's evolution scenarios, favouring 
other processes instead, such as ram pressure stripping, inter- 
nal secul ar evolution, or simple fading of the stellar popula- 
tions fsee lDresslerlll980t lAbadi et alJll999t iMoore et al.l l 1999: 
Ouilis et all 120001: lAthanassoula et ail 120051: ISmith et alJl2005t 
Tutu kov & Fedorov a 2006; Athanassoula et alj|2009alfbl) . 

Concerning simple fading, some studies support that this 
process alone cannot explain the conversion of spirals into SO's 
in clusters, as it requires the enhancement of the bulge lumi- 
nosity compared to the one o f the disk (not just its fading, see 
Christlein & Zabludoff 2004). The relevance of gas stripping in 
clusters is observa tionally proven (iKennev & Koopmannll 19991: 
IChung et al.ll2007h . but this process is not sufficient to explain 
the origin of SO's either. First, most SO's (~ 70%) reside in 
groups, not in clusters, thus ruling out this formation mecha- 
nism for the majority of them (see references below). Moreover, 
lLaurikainen et alJ J2009) argue that the stripping hypothesis (in- 
cluding quiescent star formation) cannot explain why SO's have 
a significantly smaller fraction of bars (46%) than SO/a galaxies 
(93%) or spira ls (64%-69%), which are their assume d progeni- 
tors (see also lAguerri et al.ll2009l:lBarwav et al.ll201 ll) . 

Bars themselves and internal secular processes are good can- 
didates for explaining the formation and evolution of SO's, as 
they guarantee both the enhancement of the bulge luminosity 
(inducing star formation in the centre) and a tight structural 
binding between the bulge and the disk (thanks to the inflow 
of disk material towards the galaxy centre through the bar). 
This seems to be supported by the finding of that ~ 97% of 
SO's have ovals, which might be bar relics (Laurik ainen et alJ 
2009). However, although models of internal secular evolution 
demonstrate that this mechanism can give rise to bulges of struc- 
tural properties compatible with observations of spirals, gas-free 
models (which would be the candidates to explain the formation 
of the SO bulges) produce bulges with too low concentrations 
and too large compared to real ones. Observations indicate that 
SO's have a Seisic index n ~ 2 and r e //iD ~ 0.2 on average (L10), 
whereas these model s buildup bulges with n < 1 and r e /hu ~ 0.4 
(iDebattista et al.ll2006l) . 

The environment of SO's provide clues that mergers must in- 
fluence their evolution. At z < 0.5, the fraction of SO's in groups 
is similar to that of clusters dWilman et alJl200"9|). Considering 
that most galaxies reside in groups (~ 70%, s ee iHuchra & Gelier! 
ll982tlBerlind et alJ2006t[Crook et al.ll2007l) . this means that the 
majority of SO's in the Universe are located in groups (not in 
clusters). Moreover, evolved gro ups have a much higher frac- 
tion o f SO's than young groups dFreeland et alJ 12009: Bai et al.l 
l2010l) . This morphological transformation from spirals towards 
SO's in groups is obviously driven by the processes govern- 
ing galaxy evolution in these environments, which are mergers 
and tidal interactions. Additionally, at least ~ 16% of normal 
relax ed SO's exhibit tidal debris, again suggesting past merg- 
ing dKim et alJ 120121) . and recent studies show that the low- 



level star formation activity in early-type galaxies (E-SO's) at 
intermediate redshifts is mostly driven by minor mergers too 
dKavirai et alj|2009ll201 ll) . The spread in the ages of the globu- 
lar clusters in SO's compared to ellipticals might be also point- 
ing to a relevant impact of minor mergers in their formation, 
as these processes could bring bluer (and on average younger) 
globular cl usters with them from th e smaller hosts that could 
explain it dChies-Santos et al] 1201 ll) . It is thus straightforward 
to start questioning the traditional picture of mergers as catas- 
trophic events, incapable of reproducing the structural bulge- 
disk binding observed in SO's. 

Simulations show that major mergers (i.e., encounters with 
mass ratios above 1:3) build up ellipticals, not SO-like sys- 
tems, unless the event involves a n extraordinary gas amount 
or a large bulge in the progenitors (iGonzalez-Garcfa & Balcells 
120051: iNaab & Truiillol2006llBrook et alJl2007l) . But even in this 
case, the r e //io ratios are too high (~ 0.4) compared to the typi- 
cal ones of real SO's (~ 0.2, see L10) and the remnants lie well 
below observational data in the r e -hu plane. However, interme- 
diate mergers (i.e., events with mass ratios between 1:4 and 1:7) 
can produce remnan ts that photometrically and kinema tic ally re- 
semble SO galaxies (|Be^|^^|BournMJetani^05l). 

Although simulations suggest that minor-to-intermediate 
mergers are capable of reproducing the observed distribu- 
tion of galaxies in the Seisic index (n) vs. B/D plane, very 
few studies have compared the distrib utions of the bulge 
and disk scale-lengths with observations |Aguerri et"aTI [2001: 
Scannapieco & Tissera 2003; Eliche-Moral et al. 2006a). In par- 
ticular, IScannapieco & Tisseral d2003l) studied the effects of 
mergers on the structural properties of disk-like systems in a 
cosmological context, and found that the progenitors of these 
objects at different redshifts exhibited global structures compati- 
ble with real galaxies, whereas their counterparts at z = do not 
(compare their Figs. 2 and 5). If intermediate and minor merg- 
ers have played any role in the evolution of SO's, the remnants 
of merger simulations should not just be compatible with the 
distribution of real SO's in the r s -hu plane, but they should also 
reproduce the lack of correlation found in SO's between r e /hu 
and n. The analogous result for spirals has been interpreted as 
an evidence of the fact that the Hubble Sequence is scale free 
dCourteau et alJll996l:ld7jondll996l:lGrahamll2001b1:lHunt et al.l 
2004j: lBalcells et al.ll2007l) ~ 

Therefore, we have re-addressed the question of whether in- 
termediate and minor mergers can induce evolution compatible 
with the bulge-to-disk coupling observed in real SO's or not. The 
aim of the present paper is to demonstrate that there is no reason 
to exclude dry minor mergers as a valid mechanism for driv- 
ing the evolution of lenticular galaxies along the SO sequence. 
We have analysed the evolution induced by gas-free intermedi- 
ate and minor merger events onto SO galaxies, with mass ratios 
ranging from 1:6 to 1:18. 

The paper is structured as follows. The models are briefly 
described in Sect. [2] The morphology of the final remnants is 
analysed through the simulation of realistic images in Sect. [3] In 
Sect. [4] we describe the bulge-disk decompositions performed to 
them. The structural analysis of the remnants and the comparison 
with real SO's are presented in Sect. [5] The discussion of the re- 
sults can be found in Sect. [6] Final conclusions are addressed in 
Sect. [7] A concordant cosmology is assumed for image simula- 
tions of galaxies at different cosmological distances (Qm = 0.3, 
Q A = 0.7, H = 70^8-! Mpc" 1 , see lSoergel et"aLll2007h . All 
magnitudes used in the paper are in the Vega system. 
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Table 1. Orbital and scaling parameters of each merger experiment 



Model code 


EM06 code 


M 2 /M, 


^pericenlre / 


0i n 


Primary B/D 


ar T F 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(a) M6 Ps Db 


M2TF35D 


1:6 (M6) 


0.73 (Ps) 


30 (D) 


0.5 (b) 


3.5 


(a2) M6 Ps Db TF3 


M2TF3D 


1:6 (M6) 


0.73 (Ps) 


30 (D) 


0.5 (b) 


3.0 


(a3) M6 Ps Db TF4 


M2TF4D 


1:6 (M6) 


0.73 (Ps) 


30 (D) 


0.5 (b) 


4.0 


(b) M6 Ps Rb 


M2R 


1:6 (M6) 


0.73 (Ps) 


150 (R) 


0.5 (b) 


3.5 


(c) M6 PI Db 




1:6 (M6) 


8.25 (PI) 


30 (D) 


0.5 (b) 


3.5 


(d) M6 PI Rb 




1:6 (M6) 


8.25 (PI) 


150 (R) 


0.5 (b) 


3.5 


(e) M6 Ps Ds 




1:6 (M6) 


0.87 (Ps) 


30 (D) 


0.08 (s) 


3.5 


(f) M6 Ps Rs 




1:6 (M6) 


0.87 (Ps) 


150 (R) 


0.08 (s) 


3.5 


(g) M9 Ps Db 


M3TF35D 


1:9 (M9) 


0.79 (Ps) 


30 (D) 


0.5 (b) 


3.5 


(g2) M9 Ps Db TF3 


M3TF3D 


1:9 (M9) 


0.79 (Ps) 


30 (D) 


0.5 (b) 


3.0 


(g3) M9 Ps Db TF4 


M3TF4D 


1:9 (M9) 


0.79 (Ps) 


30 (D) 


0.5 (b) 


4.0 


(h) M9 Ps Rb 


M3R 


1:9 (M9) 


0.79 (Ps) 


150 (R) 


0.5 (b) 


3.5 


(i) M18 Ps Db 


M6TF35D 


1:18 (M18) 


0.86 (Ps) 


30 (D) 


0.5 (b) 


3.5 


(j)M18PsRb 


M6R 


1:18 (M18) 


0.86 (Ps) 


150 (R) 


0.5 (b) 


3.5 


(k)M18Pl Db 




1:18 (M18) 


8.19 (PI) 


30 (D) 


0.5 (b) 


3.5 


(l)M18PlRb 




1:18 (M18) 


8.19 (PI) 


150 (R) 


0.5 (b) 


3.5 



Columns: (1) Model code: MmP[l/s][D/R][b/s][TF3/4], consult the text. The identification letters in parentheses for each model are the same as 
those used in EMI 1. The models with q;tf = 3.0 or 4.0 from EM06 have been labelled with letters "a2", "a3", "g2", and "g3" to refer them to their 
analogues with ff TF = 3.5. (2) Model code in EM06, for those models that were already presented in that paper. (3) Luminous mass ratio between 
satellite and primary galaxy. (4) First pericentre distance of the orbit, in units of the original primary disk scale-length. (5) Initial angle between the 
orbital momentum and the primary disk spin. This angle determines whether the orbit is prograde (direct) or retrograde. (6) Bulge-to-disk ratio of 
the original primary galaxy used in the experiment: B/D = 0.5 (big bulge) or B/D = 0.08 (small bulge). (7) Value of <2tf assumed for the scaling 
of the satellite to the primary galaxy. 



2. Description of the models 

We have analysed the battery of collision-les s simulations of in- 
termedi ate a nd minor mergers describe d in lEliche-Moral et all 
d2006bl) and lEliche-Moral et all d2011l) (EM06 and EMU re- 
spectively, henceforth), in which both primary and satellite are 
modelled as disk-bulge-halo galaxies with realistic density ra- 
tios (in total, 16 models). Our minor merger experiments differ 
in the considered orbital pericentre distances, the spin-orbit cou- 
pling, the primary-to-satellite luminous mass ratios, the central 
density ratios, and the galaxy types for the initial primary galaxy. 
As these models have been extensively described in the original 
papers, we just provide a brief description here, asking the inter- 
ested reader to consult the details there. 

Fourteen experiments were run using a disk galaxy with a 
prominent bulge with B/D = 0.5 as primary galaxy (equiv- 
alent to an SOb galaxy, using 185K particles in total for the 
simulation). The other two experiments have an initial primary 
galaxy with a small bulge with B/D = 0.08 (similar to an SOc, 
with 415K particles in total for the simulation). A physically- 
motivated size-mass scaling was used to ensure that the primary- 
to-satellite density ratios are realistic, forcing both galaxies to 
obey the Tully-Fisher relation: L ~ V g ^ F , being L the galax y 
luminosity and V its rotational velocity (Tullv & Fished 1 19771) . 
Different values were initially considered for the exponents of 
this relation spanning the range of observations (aiF = 3.0, 3.5, 
and 4.0). A higher ajp exponent implies a denser satellite com- 
pared to the central density of the primary galaxy. Satellites are 
scaled replicas of the primary galaxy model with a big bulge in 
all the experiments. 

The galaxy models we re built using the GalactlCS code 
dKuiiken & Dubinski|[l995l) . The disks of the primary galaxies 
and the satellites follow exponential surface density profiles ra- 
dially and vertically. Satellites and primary galaxies were al- 
lowed to relax in isolation for about ten disk dynamical times 
prior to placing them in orbit for the merger simulations, to as- 



sure that no relevant resonant structures or spiral patterns appear 
in any disks. This guarantees that our initial galaxies are SO's 
and that the evolution observed during the merger is basically 
driven by it, even the transitory bars and spirals that appear, and 
not by intrinsic disk instabilities. 

The primary galaxy with a large bulge scales to the Milky 
Way (MW) considering R = 4.5 kpc, v = 220km s _1 , and 
M = 5.1 X 1O 1O M as the units of length, velocity, and mass, 
respectively (the corresponding time unit then being 20.5 Myr). 
The primary galaxy with a small bulge matches NGC253 using 
as units R = 6.8 kpc, v = 510km s~\ and M = 2.6 x 10 n M o 
instead (in this case, the time unit is 11.7 Myr). These values, 
especially when using an appropriate M\ um /L ratio, yield total 
mass-to-light ratios close to observations in both galaxy models 
(M/L ~ 10, where M = M\ am + M dalk ). 

The galaxies are set in parabolic orbits with initial separa- 
tion equal to 15 times the primary disk scale-length in all the ex- 
periments. Two pericentre distances have been analysed: a short 
one (equal to /zq) and a long one (equal to 8/;d). The initial 
inclinations between the orbital plane and the primary galac- 
tic plane were fixed to 30° in direct orbits and to 150° in ret- 
rograde orbits. The e volution of the models was computed us- 
ing the TREECODE dHernauistlll987[ iHernauist & Katj[l989t 
lHernauisj[l^^lHernauist& Barneslll990h a nd the GADGET- 
2 codes dSpringel et aLlBOOlt ISpringelll2005l) . Using a soften- 
ing length of s — 0.02 in model units and an opening angle of 
6 = 0.6, the codes compute forces to within 1 % of those given 
by a direct summation and preserve the total energy to better 
than 0.1%. We evolved all models for ~2-4 halo crossing times 
beyond a full merger to allow the final remnants to reach a quasi- 
equilibrium state with a good conservation of energy in all runs. 
The time period between the full merger and the end of each sim- 
ulation ranges from 0.4 to 1 . 1 Gyr depending on the model, for 
the scalings commented above (see Table 2 in EM06 and Table 3 
in EMI 1 for the full-merger and total-run times of each model). 
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TableQ] lists the main characteristics of each merger experi- 
ment. Encounters with galaxy luminous mass ratios between the 
satellite and the primary galaxy equal to 1:6, 1:9, and 1:18 have 
been run. We refer to each model throughout the paper according 
to the code used in EMU: MmP[l/s][D/R][b/s], where m indi- 
cates the bulge-to-satellite mass ratio (m = 6, 9, or 18 for models 
with luminous mass ratios equal to 1:6, 1:9, 1:18, respectively), 
"PI" indicates long pericentre and "Ps" short pericentre, "D" or 
"R" describes the spin-orbit coupling ("D" for direct encounters 
and "R" for retrograde ones), and the next "b" or "s" letter indi- 
cates if the primary galaxy had a big or small bulge. We append 
to the code a final "TF3" or "TF4" suffix to include the models 
from the EM06 sample that assume q-tf = 3.0 or 4.0, respec- 
tively (excluded in EMI 1). Otherwise, a scaling with ajp = 3.5 
must be assumed. For more details, see EM06 and EMI 1. 

3. Image simulations of the remnants 

We have analysed the morphology of our remnants to check 
whether they can be considered as typical SO's or not (e.g., 
ellipticals or extremely distorted systems instead). Despite the 
definition of ver y efficient quantitative morphological indices 
dLotz et al.ll2.008l) . visual inspection is still the most trustworthy 
method to classify galaxies morphologica lly dSchawi nski et al. 
120071: Uogee et all 120091: iDarget all 1201 Ol) . So, we have clas- 
sified visually the morphology of the remnants using artificial 
images simulating the limiting magnitudes and spatial resolu- 
tions of specialized surveys. In particular, we have used the 
NIRS0S s ample as a re ference (Near-IR SO galaxy Survey, see 
lLaurikainen et al.l201 lb . As our models trace the dynamical evo- 
lution of the stellar mass, we have considered a mass-to-light ra- 
tio in the /T-band of M^^/Lk ~ 1 M b /Iq.k, which is typica l 
of early-type galaxies (Reshetnikov 2000; Portinari et al. 2004). 
Therefore, a direct equivalence between the surface density and 
the Zf-band surface brightness maps can be established in our 
models for a given mass scaling. 

We have assumed the mass and length scaling commented in 
Sect.[2] The primary galaxy in the models with big primary bulge 
is scaled to the mass and size of the MW, while the one with 
small bulge is scaled to NGC253. Note that the young stellar 
populations or dust content of these two real galaxies (which are 
spirals) are not expected to affect significantly to the stellar mass 
and typical /T-band scale-lengths of these galaxies. We have per- 
formed estimates of the /T-band magnitude of our remnants as- 
suming typical star formation histories for spiral and quiescent 
galaxies. We find that the difference in the total A'-band magni- 
tude of the galaxy between assuming quiescent or exponentially- 
decaying star formation histories typical of Sc and Sb galaxies 
is 0.06 mag at most. Therefore, if the galaxies to which we scale 
were quiescent, their luminosity and size in this band would be 
practically the same. This simple exercise guarantees a realistic 
scaling of our models in this band. 

For the MW, a radius R = 13.9kpc dMinniti et al.ll201 ll) 
and a total magnitude Kj = -23.8 mag have been adopted, 
both measured in the K band. This last value has been esti- 
mated through the convers ion of the MW ff-band tot al magni- 
tude (Bj = -20.3 mag, see iHeraudeau & Simienlfl997h to the K 
band, assuming the typical colour of nearby early-ty pe galaxies 
in the 2MASS Large Galaxy Atlas (B - K ~ 4 mag, Ijarrett et all 
12003b . For NGC253, we consider a radius R = 14.2 kpc and a 
total magnitude Kj = -23.8 mag, also in the /T-band. The ab- 
solute /f-band magnitude of this galaxy is derived assuming the 
apparent A"-band magnitude and the distance to this galaxy re- 
ported bv lJarrett etall d2003l m K = 3.77 mag at D = 3.5 Mpc). 



As indicated above, we have considered the limiting K- 
band magnitude and spatial resolution of NIRS0S in our simu- 
lated images. This sample consists of ~ 180 early-type nearby 
galaxy disks (~ 120 being SO's) with visual morphological 
classifications from the Third Reference Cata logue (RC3, see 
Ide Va ucouleur s et al ] [l9^ICorwinltZI[T99l . This survey has 
obtained multi-component decomposition of these galaxies in 
the /f-band, to ensure that the derived photometric parameters 
are not affected severely by dust and star formation. The limiting 
magnitude of the survey is ~ 21-22 mag arcsec~ 2 , and the K- 
band images have an average seeing with FWHM~ 0.7 arcsec. 
The simulated images have been convoluted by a Gaussian with 
this FWHM value to simulate the blurring effect associated to 
any PSF, and any structure below this limiting magnitude has 
been removed in the images. The typical total /f-band absolute 
magnitude of the NIRS0S galaxies is Kj ~ -24 mag, a very 
similar value to those exhibited by the two galaxies selected for 
scaling our models. 



We have considered three characteristic distances for placing 
our remnants: the minimum, average, and maximum distances of 
the galaxies in the NIRS0S sample (D ~ 6, 30, and 300 Mpc). 
The angular scale corresponding to each distance is est imated 
assuming a concordant cosmology (Spergel et all l200"7h . Both 
the cosmological dimming and the K-corrections at first order of 
approximation are included in the image simulation, albeit their 
effects are insignificant for such nearby systems. The morphol- 
ogy of the remnants at three different inclinations with respect to 
the line of sight have been analysed (face-on, inclined by ~ 45°, 
and edge-on). We have assumed that the typical noise of the K- 
band images of NIRS0S is negligible compared to the signal of 
the images relevant for morphological classification, and thus, 
no noise has been included in the images. We have checked that 
the surface brightness is conserved with distance. 

FigureQ] shows the resulting simulated images for the rem- 
nants of models M18PlDb and M6PsDs for different galaxy in- 
clinations and at the three distances considered (D = 6, 30, 
and 300 Mpc). The remnants exhibit a relaxed bulge+disk struc- 
ture with none or weak spiral arms at all distances and inclina- 
tions under consideration. This occurs in all the simulated rem- 
nants, meaning that all our dry intermediate and minor mergers 
give rise to SO-like systems, relaxed enough to be considered as 
undisturbed (the time passed since the full merger varies from 
~ 0.4 to ~ 1.1 Gyr, depending on the model, see Sect.|2|. 



No bars are observed in any models with big primary bulges, 
but those with small primary bulges present weak bar-like 
distortions and residual spiral patterns in the centre (see the 
figure). This difference arises in the lower initial B/D ratio of 
the primary galaxy, as the presence of an y central mass con- 
centratio n tends to stabilize t he di sk dPfenniger & Norman! 
Sellwood& Moord Il999t 
Athanassoula et al. 2005; 
Gonzalez-Garcfa & Balcelisl 120051: ICoxetal.1 l2008h . This 
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2002 



Bournaud & Combes 
Deb attista et al] 12005: 



makes spiral arms and bar distortions to be more efficiently 
inhibited in the remnants of models with initial big primary 
bulges than in those with small ones. Nevertheless, these inner 
structures are insignificant compared to the bulge and disk 
component, so we can assume that the luminous mass in our 
remnants is distributed basically according to an axi-symmetric 
bulge+disk structure. 
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Fig. 1. Simulated Zf-band images of the final remnants of models M18PlDb (upper figure, experiment k in Table[T} and M6PsDs (lower figure, experiment e in the same Table). 
We show different galaxy views for the two further simulated distances. The largest frame on the left of each figure corresponds to a face-on view of the remnant at D = 6 Mpc. 
Surface brightness units in the K band are in mag arcsec -2 . A limiting magnitude of = 22 mag arcsec -2 and a spatial resolution of FWHM= 0.7 arcsec have been assumed 
(see the text for more details). A grey colour palette is used to represent the different surface brightness levels, set according to a logarithmic scale to enhance spiral patterns 
in the images with respect to the central bulge. Whiter colours indicate brighter regions (consult the levels bars). The size of each frame fits exactly the apparent angular size 
of the galaxy for the distance considered in each case. There are yellow squares at the top left corner of each frame indicating the size of the resolution element in each image 
(0.7x0.7 arcsec 2 ). 
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Table 2. Bulge Seisic indices, bulge-to-disk luminosity, and scale-length ratios of the original primary galaxies and final remnants 
of the merger experiments 



Model 




n 


B/D 


r e /h D 


(1) 




(2) 


(3) 


(4) 


Original galaxy, B/D = 


0.5 


0.92±0.12 


0.51+0.03 


0.190+0.006 


Original galaxy, B/D — 


0.08 


0.52±0.02 


0.0760±0.0024 


0.1611+0.0005 


(a) Mo rs Do 




1. /+U.4 


0.8Z+0.1U 


U.164+U.UU5 


(azj Mo rs Do Irj 




i.o+u.i 


n hh i n no 
U. / /+U.UB 


U. 1 /+0.01 


(a j ) Mo r S DO 1 r4 




1 O i Ci £ 
1.8+1). 


U.9+U.3 


U. 109+U.U1 1 


(b) M6 Ps Rb 




1.7±0.4 


0.70+0.10 


0.182+0.008 


(c) M6 PI Db 




1.11 ±0.04 


0.676+0.024 


0.2396+0.0002 


(d) M6 PI Rb 




1.02±0.06 


0.55+0.03 


0.2734+0.0014 


ye ) IVlO rS US 




Z.ZV+U. 10 


n 1 ^Q-i-n m i 
U, 1 /y+u.uiz 


U. lDU±U.UUo 


(f) M6 Ps Rs 




2 44+0 22 


140+0 012 


143+0 010 


(g) M9 Ps Db 




1.5+0.3 


0.770+0.08 


0.160+0.005 


(g2) M9 Ps Db TF3 




1.3+0.3 


0.65+0.06 


0.166+0.009 


(g3) M9 Ps Db TF4 




1.3+0.3 


0.75+0.08 


0.174+0.004 


(h) M9 Ps Rb 




1.9±0.8 


0.72+0.18 


0.165+0.015 


(i)M18Ps Db 




1.5+0.3 


0.60+0.08 


0.157+0.013 


(j)M18PsRb 




1.6±0.3 


0.68+0.07 


0.159+0.006 


(k)M18Pl Db 




1.05±0.05 


0.60+0.03 


0.2206+0.0002 


(l)M18PlRb 




1.06±0.05 


0.51+0.03 


0.1974+0.0004 



Columns: (1) Model code from TablefT] (2) bulge Sersic index, (3) bulge-to-disk mass ratio derived from the Sersic+exponential fit (equivalent to 
bulge-to-disk luminosity ratio in the A*-band, see Sect.f3](, (4) bulge-to-disk scale-lengths ratio (defined as the ratio of the bulge effective radius to 
the disk exponential scale-length). 

whole radial range, a quite reasonable result compared to typical 
observational errors. Errors of the p hotometric parameters are 
estimated w ith the bootstrap method (lEfron & Tibshiranil [l993: 
iPress et al]|1994t) . For all our experiments, we recover a two 
component system with the outer parts well fitted by an expo- 
nential profile (i.e., a disk), while the inner parts are well repre- 
sented by Sersic functions with 1 < n < 2.5 (see Fig. 3 in EM06 
for some examples). Final values of n, B/D, and r e //io derived 
from these fits to the surface density profiles of the remnants are 
listed in Table [2] Our models are scalable in size, because they 
are collision-less. This means that the bulge effective radius of 
each remnant can be scaled to any physical size, just considering 
that the disk scale-length is related to this value according to the 
r^/hu ratio of this remnant tabulated in Table|2] 

5. Results 

We have analysed the evolution induced by the mergers in the 
bulge and disk photometric parameters of the primary galaxy, 
to compare the final values of B/D, r e /hu, and n in the rem- 
nants with the observational distri butions of these par ameters in 
real spirals and SO's reported bv lWeinzirl et alj ([2009) and L10. 
These studies use images in NIR bands (H and K, respectively). 
We also compare our results with the evolution entailed by ac- 
cretion experiments onto an SOb of high-density dSph satellites 
(Agu erri et al.ll200 l |) and with internal secular models of pure 
exponential disks (iDebattista etlil1l2006l) . 



5. 1. Evolution along the sequence of SO Hubble types 

The initial B/D ratios of the primary galaxies increase after 
the merger by a factor ~ 1.1-1.8 in the case of big primary 
bulges and ~ 1.8-2.3 in the models with small ones (see Table|2]i. 
According to the observe d correlations between n and B/D wit h 
the morphological type dGrahamll2001bt IWeinzirl etakl 12009). 
this means that the remnants of the experiments starting with an 



4. Bulge-disk photometric decompositions of the 
remnants 

We have derived the bulge-disk photometric parameters from 
azimuthally-averaged surface brightness profiles for the lumi- 
nous matter, accounting for the axial symmetry of all the rem- 
nants (see Sect. [3]). In general, the global photometric param- 
eters derived from ID and 2D fits in axi-symmetric case s are 
quite similar within errors (Ide Jong|ll996l lMahaballll998l) . The 
ID profiles present a much better signal-to-noise as compared 
to the 2D s urface brightness maps, ensuring a higher accuracy 
of the fits dPenpetal.1 12010). We have derived the ID surface 
brightness profiles with all the merger remnants viewed face-on 
to compare with observations, which usually correct their results 
for the effects of galaxy inclination. 

We have fitted a combined Sersic+exponential function to 
the density profiles, as in EM06. The Sersic profile has been ex- 
tensively use d for modell ing bulge and elliptical surface bright- 
ness profiles (ISersidll968l) : 

I(r) = I e exp{-b n [(r/r e ) l/ "-l]}, (1) 

where r e is the bulge effective radius, I e is the surface brightness 
at r e , and n is the Sersic index, which is an indicator of the bulge 
concentration. The factor b n is a function of n, which may be 
approximated by b„ = 1 .9992 n -0.3271 in t he range 1 < n< 10, 
with an error < 0.15% in the whole range (ICiottilll991tlGrahaml 
[2001ah . 

The expon ential law adeq uately describes the radial profiles 
of most disks (lFreemanll9 70): 

/(r)=7 exp(-r//! D ), (2) 

/id being the scale-length of the profile and Io representing the 
extrapolated central surface brightness. 

The code used for fitting the radial surface density profiles 
is described in EM06. Residuals are below ~0.2 mag along the 
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Fig. 2. Growth vectors in the log («)-log (B/T) plane driven by 
our merger experiments, compared to the observational distribu- 
tions of SO's and spirals. Each segment starts at the location of 
the original galaxy model and ends at the final n and B/T val- 
ues of the remnants. The growth vectors corresp onding to the 
merge r models with dense spheroidal satellites bv lAguerri et al.l 
(1200 ll) are represented too, as well as the bulges resulting from 
th e internal secular evolu tion models in pure exponential disks 
bv lDebattista et al.l (120061) . Consult the legend in the figure. 



SOb primary galaxy have become SOa's, while the SOc primary 
galaxies are transformed into SOb's. 

The dry intermediate and minor mergers lead to a more no- 
ticeable growth of the bulge-to-total luminosity ratio (B/T) in 
the experiments with small primary bulges (~ 150% with re- 
spect to the original value) than in those with big primary bulges 
(~ 60%), for similar mass ratios and orbits of the encounter (con- 
sult Table©. 

Figure© shows the final distributions of all remnants in the 
log (n)-\og(B /T) plane, compared to the observed distributions 
of real galaxies with different morphological types. For each 
merger simulation, we plot the initial and post-merger values. 
The figure shows that dry intermediate and minor mergers onto 
SO's increase the mass (B/T) and concen tration of the bulge 
(n), as already repor ted by previous studies ( Ag uerri et al.l 2001: 
Younger et al. 2007). Although the initial SO models are offseted 
from the location of typical SO's in this plane (the original n 
values are lower than it should be for the B/T ratio), the rem- 
nants are located in the region populated by SO's. The evolu- 
tion induced by these mergers along the sequence of SO Hubble 
types is coherent with the observational fact that the majority 
of SO's have exponential-like bulges (1 < n < 2), instead of 
n > 3 bulges (lLaurikainen et al.l2005. 2006; Balc ells et all2 007: 



Wein zirl et al.l l2009). The bulges resulting from i nternal secular 
evolution as modelled by Debatti sta et al.l (120061) overlap in the 
plane with our remnants. 
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Fig. 3. Growth vectors in the log (hY))-\og (r e ) plane driven by our 
merger experiments, compared to the observed distributions of 
SO's and spirals. The straight lines show the correlations found 
by L10 for different morphological types (red solid line: for the 
SO's, blue dotted-dashed line: for the S0/a+Sa's, green dashed 
line: for the Sab-Sbc's). The blue shaded region indicates the lo- 
cations in the plane that are covered by our models adjusting the 
scaling (see Sect.©. Note that the change in both scale-lengths 
is negligible, making the remnants to basically overlap with the 
original location of the primary galaxies in the plane. Consult 
the legend for the observational data in the figure. The legend 
for the models is the same as in Fig.© 

Figure© shows some trends of the bulge growth driven by 
the merger with the initial conditions of the model: 

1 . The accretion of a satellite by a galaxy with an original big 
bulge leads to a smoother growth in both n and B/T than in 
the case of being accreted by a galaxy with a small bulge, 
relatively to their initial values. 

2. Moreover, the minor mergers undergone in orbits with long 
pericentres produce a much more moderate growth in both n 
and B/T than analogous mergers with short impact parame- 
ters. 

3. Curiously, the evolution in the n-B/T plane seems to depend 
more significantly on the orbit than on the mass ratio of the 
encounter for the same original gala xy, contrary to previou s 
findings in major merger simulations (iBournaud et al.l2004l) . 

4. Finally, as already reported in EM06, mergers with denser 
satellites (i.e., higher o-tf) lead to more noticeable growths 
in « than those with low-density sate llites. The high-density 
satellites of the lAguerri et al.l d2001l) models produce more 
extreme evolution in this plane than the one produced by 
our models for the same original galaxy. We expect satel- 
lites with lower central densities than ours to disrupt com- 
pletely before reaching the galaxy and contribute to the stel- 
lar ha lo of the galaxy, as already observed in other simula- 
tions dCooperetal.ll2010h . 

Our models suggest that dry intermediate and minor mergers 
can induce bulge growth in two ways relatively to their initial 
state, depending on the mass of the original primary bulge. On 
the one hand, galaxies with small bulges experience a promi- 
nent rise in «, and a much more moderate increase of B/T. This 
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is because the satellite core ends deposited in the remnant cen- 
tre without disrupting, rising the mass concentration, and thus n 
(EMI 1). On the other hand, galaxies with big bulges experience 
moderate growths in n and B/T after the accretion of a satellite 
compared to their original status, with a higher dependence on 
the orbit than on the mass ratio. This is because satellites disrupt 
completely before reaching the centre, depositing their mass at 
longer distances from it (EM06; EMI 1). 

In conclusion, these models confirm that dry minor and in- 
termediate mergers onto SO's can induce evolution in the n-B/T 
plane towards earlier SO Hubble types. 

5.2. Correlations between the bulge and disk scale-lengths 

In Fig. |3] we show the location of the final remnants of our 
merger models in the log (/!o)-log (r e ) plane, compared to the 
observational distributions of real galaxies with different mor- 
phological types. We have marked the locations in the plane that 
would be achievable by our models just assuming a different 
scaling. The straight lines indicate the correlations in the plane 
found for different galaxy types by L10. Observationally, galax- 
ies distribute diagonally in this diagram, i.e., the scale-lengths of 
the bulge and disk components correlate tightly in real galaxies. 
These correlations are quite similar for different galaxy types 
(compare the straight lines) and, more relevantly, the distribu- 
tions of the different types overlap completely in the diagram. 
The figure shows that the region in the h-Q-r^ plane that can be 
achieved by our models through different scaling overlaps with 
the distribution of real SO's basically. 

Our models demonstrate that dry intermediate and minor 
mergers can produce a negligible change in the location of a 
galaxy in the hu-r e plane, independently on the B/T of the orig- 
inal SO and on the properties of the encounter. Therefore, this 
sort of mergers basically contributes to rise slightly the disper- 
sion of the relation between the scale-lengths of the bulge and 
disk components of SO's, although they can induce a noticeable 
bulge growth and move the galaxies towards earlier types. This 
result indicates that dry intermediate and minor mergers are ca- 
pable of producing remnants in which the scale-lengths of the 
bulge and the disk are coupled, if they were already coupled in 
the original galaxy. Hence, the main conclusion of Fig.[3]is that 
minor mergers cannot be discarded as feasible mechanisms for 
the evolution of SO's just on the basis of the strong coupling of 
the bulge and disk scale-lengths observed in these galaxies. 

We plot in Fig.|4]the evolution of r e and hu versus n induced 
by our minor merger experiments compared to the observational 
distributions of different morphological types. The blue shaded 
region indicates the locations of these diagrams that are achiev- 
able with our models, just assuming a length scaling different to 
the one indicated in Sect.[2] The simulations indicate that the dis- 
tribution in the plane of SO's could be reproduced through dry in- 
termediate and minor merger events basically, starting from SO's 
of a late Hubble type. Therefore, Fig s.[2]|4] prove that the global 
bulge-to-disk structural properties of real SO's are not incompat- 
ible with an evolution driven by gas-poor mergers of low mass 
ratio. 

Figure|5] shows the evolution of r e and versus B/T in- 
duced by our minor merger experiments compared to the ob- 
servational distributions of different morphological types. Both 
scale-lengths do not vary by more than ~ 20% of their origi- 
nal value after the merger. The final B/T ratio tends to rise after 
the mergers, except in a few cases bv lAguerri et al.l (|2001). The 
simulations can reproduce the distribution of SO's also in these 
planes, just assuming a different scaling. 



In conclusion, our models demonstrate that the dry interme- 
diate and minor mergers can evolve SO's along the sequence of 
SO Hubble types towards earlier types, fulfilling global bulge-to- 
disk structural relations compatible with observations. 

5.3. A scale-free Sequence of SO Hubble types 

The lack of correlation of r e /hu with n and B/T for spirals 
led many authors to conclude that the Hubble sequence is 
scale-free dde Jong|ll996t ICourteau et al.lll996t; lGraharnll2001al: 
iMacArthur et alJl2003l: iBalcells et al.ll2007h ."Figurel6lrepresents 
the evolution driven by our minor merger experiments in r e /hu 
as a function of the Sersic index (left panel) and B/T ratio (right 
panel), compared to observational data. The data support the pre- 
vious finding also for SO's: there is a wide dispersion in the rjh-o 
ratio for any n and B/T values found in SO's too, although this 
dispersion rises at the region of high n and B/T values, where 
SO's tend to accumulate (this is more noticeable in the right 
panel). 

The remnants of our experiments are also compatible with 
the observed distributions of SO's in these planes. Our models 
show that dry mergers of low mass ratios onto SO's make them 
to move towards higher n and B/T values, keeping the origi- 
nal rjh-o ratio. Although the minor merger can produce a sig- 
nificant rise in n or B/T depending on the initial conditions of 
the encounter, the ratio r t /h-Q is preserved or slightly dispersed. 
This demonstrates that the evolution induced by these mergers is 
compatible with three observational facts: 1) with the accumula- 
tion of SO's of early-types towards high n and B/T; 2) with the 
lack of correlation of r e /ho with the morphological type of the 
SO, because this ratio depends exclusively on the original galaxy 
experiencing the accretion, not on the minor merger itself; and 
3) the slight increase of dispersion in r e /hu as the Hubble type 
becomes earlier. 

Note that the i nternal secular evolutio n models of pure ex- 
ponential disks by Debattist a et al.l (120061) generate bulges that 
do not overlap with real SO's in Fig. [6] even those without gas. 
Their structures are more similar to the bulges of spirals than to 
those of SO's. 

Therefore, our models suggest that dry mergers of low mass 
ratios onto SO's can provide one of the feasible explanations for 
the buildup of a scale-free sequence of Hubble SO types, starting 
from SOc's of late types formed through other mechanisms. 

6. Discussion 

Our simulations prove that dry intermediate and minor mergers 
can make SO galaxies to evolve towards earlier SO types, ful- 
filling global structural relations observed between the bulges 
and disks of these galaxies. According to the final B/T ratios of 
the remnant SO's (see Tabled, our experiments show that these 
mergers induce the following transformations: SOc — >S0b and 
SOb — >S0a. These evolutionary paths are indicated with green 
arrows in the Hubble Sequence diagram of Fig. [7] 

Considering that the growth in n and B/D of successive 
minor mergers is cumulative (although it seems to saturate at 
n ~ 4-5 values, see lAguerri et al.l l2001; Bou rnaud et al.ll2007h . it 
should be possible to reproduce the distributions of SO's in the 
structural planes shown in Figs. |2]|6] through subsequent minor 
merging. Therefore, the models demonstrate that dry intermedi- 
ate and minor mergers can give rise to a sequence of SO Hubble 
types parallel to the one of spirals, starting from SOc's. 

As commented in Sect.Q] mergers are usually discarded as 
feasible evolutionary mechanisms for SO's, attending to several 
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Fig. 4. Growth vectors in the log (r e )-log («) and log (ZioHog (n) planes driven by our merger experiments, compared to the obser- 
vational distributions of SO's and spirals. The blue shaded region indicates the locations in the plane that are covered by our models 
just using a scaling different to the one indicated in Sect. [2] Consult the legend for the observational data in the figure. The legend 
for the models is the same as in Fig. [2] 
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observational distributions of SO's and spirals. The blue shaded region indicates the locations in the plane that are covered by our 
models just using a scaling different to the one indicated in Sect. [2] Consult the legend for the observational data in the figure. The 
legend for the models is the same as in Fig. [2] 



observational facts that people interpret as exclusive of pure in- 
ternal secular evolution. By the contrary, the present models 
prove this interpretation to be questionable, as many of these 
observational results can be explained or are compatible with an 
evolution driven by dry intermediate and minor mergers. 

In particular, the fact that the bulge and disk scale-lengths in 
SO's are tightly correlated is compatible with an evolution driven 
through these mergers, because they nearly preserve the origi- 
nal bulge-to-disk scale-lengths ratio. This evolutionary scenario 
only requires the correlation between the two scale-lengths to 
have been established at an early epoch of bulge growth. But 
this is ensured by observational data, as L10 demonstrate that 
this correlation does not depend on galaxy type. This means 
that it is completely feasible to partially build up the parallel 
sequence of SO Hubble types through dry minor mergers onto 



SO's, as these processes modify the sizes of both components in 
such a way that their ratio is nearly conserved, moving galaxies 
towards earlier Hubble types at the same time. Moreover, these 
mergers introduce some dispersion in the r e //io ratio compatible 
with observational data, providing a straightforward explanation 
for the lack of trend of the r e /ho with the Hubble type in SO 
systems. 

Curiously, another observational result compatible with a 
dry minor merger evolution is the high fraction of pse udob- 
ulges found in SO's (lLaurikainen et al.ll2005l I2006L l2007h . Note 
that the SO's resulting from our dry minor mergers have n < 2 
and B/T < 0.4. Although the definition of pseudobulges take 
into account several aspects of the shapes and dynamics of the 
central regions of galaxies, the bulges obtained in our simula- 
tions fulfill some criteria that define a pseudobulge, according 
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Fig. 7. Schematic representation of the evolutionary paths driven by our dry intermediate and minor merger experiments along the 
sequence of SO Hubble types. Bulge-to-total luminosity ratios increase towards the left, accordingly to the evolution produced by 
our simulated mergers. Our models show that minor mergers can induce the evolution of SOb — >S0a and of SOc — >S0b, reproducing 
the globa l bulge-to-disk structural re lations found in real SO's. The revised Hubble tuning fork plotted in the figure is inspired in 
Fig. 1 bv lKormendv & Bender! d2012l) . 



to iKormendv & Kennicuttl ((2004). Therefore, dry minor merg- 
ers cannot be discarded as possible evolutionary mechanisms 
of SO's just arguing that many of them exhibit pseudobulges, 
assuming that these bulge structures are associated exclusively 
with pure internal secular evolution. Our models prove that dry 
minor mergers can generate pseudobulges if the original SO had 
already one. 

In fact, our minor merger experiments show that these events 
can induce internal secular evolution, even in the absence of gas 
and star formation (EM06; EMU). In EM06, we reported that 
the bulge growth in these experiments was basically due to the 
injection of disk material towards the galaxy centre through the 
transient disk distortions induced by the interaction. The satel- 



lite material deposition was essential for rebuilding an exponen- 
tial disk at central and intermediate radii in the remnant, but not 
for the bulge growth. These transitory distortions (such as spi- 
rals, bars, and ovals) are usually ascribed to pure internal secular 
evolution, but these models prove that minor and intermediate 
mergers can be highly efficient at inducing them, even in gas- 
poor events. 

Additionally, the lower frequency of bars detected in SO's 
compared to the one found in Sa+SO/a's is not well understood 
in an evolutionary sc enario majoritarily driven by pure inter- 
nal secular processes dKormendvll979tlLaurikainen et al.f 2007 , 
[20091 iButa etalJ 120101 iNair & AbrahaBHoiOD Sa+S0/a's are 
the most probable progenitors of SO's according to this sce- 
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nario, and thus, their bar fractions should be similar, unless 
the process entails a bulge growth enough to inhibit bar distor- 
tions. But this does not seem to be the case. According to esti- 
mates by Kormendv & Kennicutt (2004), internal secular evolu- 
tion can give rise to bulges with B/T ~ .16 at most, whereas 
SO's have B/T ~ 0.25 - 0.30 typically dBalcells etal] 120071: 
lLaurikainen et al.ll2007t IWeinzirl et al.ll2009l) . Nevertheless, the 
difference in bar fractions between these two galaxy types is co- 
herent with an SO evolution partially driven by dry minor merg- 
ers. We have shown that mergers are quite efficient rising the 
concentration and mass of the central bulge, even for gas-free 
encounters with mass ratios as low as 1:18. This process tends 
to suppress bar formation, as commented above, being compat- 
ible with a decrease of the bar fraction in SO's compared to any 
other Hubble type. 

The high fraction of dynamically-cold inner components 
found in SO's is another observational result that is usually 
put forward against the possible merger-driven evolution i n 
these galaxies (see references in iKormendv & Bended 120121) . 
Although we have not focussed the present study on this, it is 
remarkable that our merger models also show this concept to be 
questionable. In EMI 1, we show that all our remnants host thin 
rotationally-supported inner components made out of disrupted 
satellite material, such as inner disks, inner rings, pseudo-rings, 
nested inner disks, and central spiral patterns. Consequently, the 
high number of SO's hosting dynamically-cold stellar compo- 
nents might be pointing to their merger-related origin instead of 
to the contrary fact, as usually interpreted. 

Moreover, the high number of ovals and lenses found in SO's 
(~ 97%) is considered as a result pointing to the relevance of 
bars (and thus, of pur e internal secular evoluti on) in the forma- 
tion of these systems (lLaurikainen et al.ll2009l) . But, again, this 
result does not exclude a minor-merger origin for these struc- 
tures. Many of our remnants have ovals, probably derived from 
the transitory bar-like distortions induced by the encounter in the 
original galaxy disk, as we also show in EMI 1. 

Finally, the apparently non-distorted morphology exhibited 
by an SO cannot be considered as a sign excluding recent past 
merging either; first, because very deep imaging may reveal faint 
tidal featu res that are not detectable in typi cal images in many 
cases (see iDuc et al. 1 1201 U iKim et al .112012). and secondly, be- 
cause this study proves that remnants of dry intermediate and 
minor mergers can exhibit a relaxed morphology only l-2Gyr 
after the complete merger for the typical observing conditions of 
current surveys. 

We must remark that all these results do not exclude the rel- 
evant contribution of other processes to the buildup of SO's, but 
the present study clearly rehabilitate dry minor mergers as feasi- 
ble mechanisms for explaining partially their evolution. 



7. Conclusions 

Our collision-less N-body simulations show that intermediate 
and minor mergers are suitable mechanisms t o drive bulge- 
growth along the SO sequence ea rly reported by van d en Bereh 
(1 19761) and recently updated by IKormendv & Bended (120121) . 
Different orbits, density ratios, and mass ratios (1:6, 1:9, 1:18) 
are considered, as well as two different models for the primary 
galaxy (SOb and SOc). 

We show that all remnants would be classified as undis- 
turbed SO's by current specialized surveys. The merger events 
induce noticeable bulge growth in all cases, inducing the follow- 
ing morphological transformations: SOc — >S0b and SOb — >S0a. 
Depending on the initial B/D of the primary galaxy that accretes 



the satellite, the merger induces bulge growth in two ways. In the 
case that the primary galaxy has a small bulge, the minor merger 
produces a prominent growth in n and a noticeable rise in B/T 
compared to their initial values, whereas in the case that it has 
a massive one, the increases in n and B/T are much more mod- 
erate for the same satellite. We also find that mergers with long 
pericentre orbits produce a much more moderate bulge growth 
than in analog cases with short impact parameters. In our exper- 
iments, the evolution in the n-B/T plane depends more signifi- 
cantly on the orbit than on the mass ratio of the encounter for the 
same original SO galaxy. 

In all cases, the merger induces a negligible change of r e and 
hu in the galaxy, in such a way that the ratio of these two scale- 
lengths is nearly preserved. Therefore, if there is a coupling be- 
tween the scale-lengths of the bulge and disk components prior 
to the merger, the merger does not break this coupling, but just 
disperses it slightly. This fact provides a simple explanation for 
the lack of correlation between the ratio of the bulge and disk 
scale-lengths and the SO Hubble type reported by observations 
for SO's. 

Our models prove that dry intermediate and minor mergers 
onto SO galaxies can make them to evolve within the sequence of 
SO Hubble types towards earlier types, fulfilling global bulge-to- 
disk structural relations compatible with those observed in real 
SO's. This means that minor mergers should not be discarded 
from the evolution scenarios of SO's just on the basis of the 
strong correlations observed between the bulge and disk scale- 
lengths of these galaxies. 
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